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ABSTRACT 

It is shown how to decorrelate the (prewhitened) power spectrum measured from a 
galaxy survey into a set of high resolution uncorrelated band-powers. The treatment 
includes nonlinearity, but not redshift distortions. Amongst the infinitely many pos- 
sible decorrelation matrices, the square root of the Fisher matrix, or a scaled version 
thereof, offers a particularly good choice, in the sense that the band-power windows 
are narrow, approximately symmetric, and well-behaved in the presence of noise. We 
use this method to compute band-power windows for, and the information content 
of, the Sloan Digital Sky Survey, the Las Campanas Redshift Survey, and the IRAS 
1.2 Jy Survey. 

Key words: cosmology: theory - large-scale structure of Universe 



1 INTRODUCTION 



Accurate measurements of the principal cosmological pa- 
rameters now appear to be within reach (Turner 1999). 
Large redshift surveys of galaxies, notably the Two-Degree 
Field Survey (2dF) (Colless 1998; Folkes et al. 1999) and the 
Sloan Digital Sky Survey (SDSS) (Gunn & Weinberg 1995; 
Margon 1998), should be a gold mine of cosmological infor- 
mation over wavenumbers k ~ 0.003-300 /iMpc -1 . Informa- 
tion from such surveys on large scales (k ^ 0.1/iMpc - ) 
should improve greatly the accuracy attainable from up- 
coming Cosmic Microwave Background (CMB) experiments 
alone (Eisenstein et al. 1999), but in fact most of the infor- 
mation in these surveys is on small scales, in the nonlinear 
regime (Tegmark 1997b). 

While the power spectrum may be the ideal carrier of 
information at the largest, linear scales where density fluc- 
tuations may well be Gaussian, it proves less satisfactory 
at moderate and smaller scales, since nonlinear evolution 
induces a broad covariance between estimates of power at 
different wavenumbers, as emphasized by Meiksin & White 
(1999) and Scoccimarro, Zaldarriaga & Hui (1999). Ulti- 
mately, one can imagine that there exists some kind of map- 
ping that translates each independent piece of information in 
the (Gaussian) linear power spectrum into some correspond- 
ing independent piece of information in the nonlinear power 
spectrum, or some related quantity. That such a mapping 
exists at least at some level is evidenced by the success of 
the analytic linear <-> nonlinear mapping formulae of Hamil- 
ton et al. (1991, hereafter HKLM) and Peacock & Dodds 
(1994, 1996). However, if such a mapping really existed, that 
was not only invertible but also, as in the HKLM-Peacock- 



Dodds formalism, mapped delta-functions of linear power 
at each linear wavenumber into delta-functions of nonlinear 
power (or some transformation thereof) at some possibly dif- 
ferent nonlinear wavenumber, then that mapping ought to 
translate uncorrelated quantities in the linear regime - pow- 
ers at different wavenumbers - into uncorrelated quantities 
in the nonlinear regime. The fact that the nonlinear power 
spectrum is broadly correlated over different wavenumbers 
shows that the HKLM-Peacock-Dodds formalism cannot be 
entirely correct. 

In the preceding paper (Hamilton 2000, hereafter Pa- 
per 3), it was shown that prewhitening the nonlinear power 
spectrum - transforming the power spectrum in such a way 
that the noise covariance becomes proportional to the unit 
matrix - substantially narrows the covariance of power. 
Moreover, this narrowing of the covariance of power occurs 
for all power spectra tested, including both realistic power 
spectra and power law spectra over the full range of indices 
permitted by the hierarchical model. It should be empha- 
sized that these conclusions are premised on the hierarchi- 
cal model for the higher order correlations, with constant 
hierarchical amplitudes, and need to be tested with iV-body 
simulations. 

In the meantime, if indeed there exists an invertible 
linear «-» nonlinear mapping of cosmological power spectra, 
then it would appear that the prewhitened nonlinear power 
spectrum should offer a closer approximation to the right 
hand side of this mapping than does the nonlinear power 
spectrum itself. Whatever the case, the prewhitened non- 
linear power spectrum has the practical benefit of enabling 
the agenda of the present paper - decorrelating the galaxy 
power spectrum - to succeed over the full range of linear 
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to nonlinear wavenumbers. That is to say, if one attempted 
to decorrelate the nonlinear power spectrum itself into a 
set of uncorrelated band-powers, then the band-power win- 
dows would be so broad, with almost canceling positive and 
negative parts, that it would be hard to interpret the band- 
powers as representing the power spectrum in any mean- 
ingful way. By contrast, the covariance of the prewhitened 
nonlinear power is already narrow enough that decorrelation 
into band-powers works quite satisfactorily. 

Paper 3 showed how to construct a near approximation 
to the minimum variance estimator and Fisher information 
matrix of the prewhitened nonlinear power spectrum in the 
Feldman, Kaiser & Peacock (1994, hereafter FKP) approxi- 
mation, valid at wavelengths short compared to the scale of 
the survey. In the present paper we describe how to complete 
the processing of the prewhitened nonlinear power spectrum 
into a set of decorrelated band-powers that come close to 
fulfilling the ideals of (a) being uncorrelated, (b) having the 
highest possible resolution, and (c) having the smallest pos- 
sible error bars. 

The idea of decorrelating the power spectrum was pro- 
posed by Hamilton (1997; hereafter Paper 2), and was fur- 
ther discussed and successfully applied to the CMB by 
Tegmark & Hamilton (1998). 

Like Paper 3, the present paper ignores redshift distor- 
tions and other complicating factors, such as light-to-mass 
bias, on the Alfa-Romeo (1969 P159 Shop Manual) principle 
that it is best to adjust one thing at a time. 



2 DATA 

This paper presents examples mainly for two redshift sur- 
veys of galaxies, the Sloan Digital Sky Survey (SDSS; Gunn 
& Weinberg 1995; Margon 1998) and the Las Campanas 
Redshift Survey (LCRS; Shectman et al. 1996; Lin et al. 
1996). The information content of the prewhitened power 
spectrum of the IRAS 1.2 Jy survey (Fisher et al. 1995) is 
also shown for comparison in Figure m 

Calculating the Fisher matrix of the prewhitened power 
spectrum, as described in §6 of Paper 3, involves evaluating 
the FKP-weighted pair integrals R(r; /i) (Paper 3, eq. [84]), 
commonly denoted {RR), of the survey. For SDSS and IRAS 
1.2 Jy, the pair integrals R(r; fj,), were computed in the man- 
ner described by Hamilton (1993), over a logarithmic grid 
spaced at 0.02 dex in pair separation r and 0.1 dex in FKP 
constant /x. For LCRS, we thank Huan Lin (1996, private 
communication) for providing FKP-weighted pair integrals 
R(r; jj,) on a grid of separations r and FKP constants fi. 

The geometry of SDSS is described by Knapp, Lupton 
& Strauss (1997, 'Survey Strategy' link). We consider only 
the main galaxy sample, not the Bright Red Galaxy sample. 
We thank D. H. Weinberg (1997, private communication) 
for providing the anticipated radial selection function of the 
survey, and M. A. Strauss (1997, private communication) 
for a computer file and detailed explanation of the angular 
mask. The angular mask of the northern survey consists of 
45 overlapping stripes each 2? 5 wide in declination, forming 
a roughly elliptical region extending 110° in declination by 
130° in the orthogonal direction. The angular mask of the 
southern survey consists of 3 separate 2? 5 wide stripes cen- 



tred at declinations 15°, 0°, and —10°. The survey covers a 
total area, north plus south, of 3.357 sr. 

The LCRS (Shectman et al. 1996) covers six narrow 
slices each approximately 1?5 in declination by 80° in right 
ascension, three each in the north and south galactic caps. 

The angular mask of the IRAS 1.2 Jy survey (Fisher et 
al. 1995) is the same as that of the 2 Jy survey (Strauss et 
al. 1990, 1992), and consists of the sky at galactic latitudes 
\b\ > 5°, less 1465 approximately 1° x 1° squares, an area 
of 11.0258 sr altogether. We thank M. A. Strauss (1992, pri- 
vate communication) for a computer listing of the excluded 
squares. The radial selection function was computed by the 
method of Turner (1979). Distances were computed from 
redshifts assuming a fiat matter-dominated cosmology, and 
galaxy luminosities L were corrected for evolution with red- 
shift 2 according to L oc 1 + z. 



3 DECORRELATION 
3.1 Decorrelation matrices 

Let F be the Fisher information matrix (see Tegmark, Tay- 
lor & Heavens 1997 for a review) of a set of estimators 6 
of parameters 8 to be measured from observations. Below 
we will specialize to the case where the parameters are the 
prewhitened power spectrum, but for the moment the pa- 
rameters 6 could be anything. Assume, thanks to the central 
limit theorem or otherwise, that the covariance matrix of the 
estimators 9 is adequately approximated by the inverse of 
the Fisher matrix 

(A§A6 T ) = F' 1 . (1) 

A decorrelation matrix W is any real square matrix, not 
necessarily orthogonal, satisfying 

F = W T AW (2) 

where A is diagonal. The quantities \N9 are uncorrelated 
because their covariance matrix is diagonal 

W(A(9A<9 T )W T = A" 1 . (3) 

Being inverse variances of W#, the diagonal elements of the 
diagonal matrix A are necessarily positive. The Fisher ma- 
trix of the decorrelated quantities \N0 is the diagonal matrix 
A, since 

(W(A0A0 T )W T )" 1 = A . (4) 

Without loss of generality, the decorrelated quantities 
W# can be scaled to unit variance by multiplying them by 
the square root of the corresponding diagonal element of A. 
Scaled to unit variance, the decorrelation matrices satisfy 

F = W T W . (5) 

There are infinitely many distinct decorrelation matrices 
satisfying equation (g). Any orthogonal rotation of a 
decorrelation matrix W yields another decorrelation matrix 
V = OW, since V T V = W T T OW = W T W = F. Con- 
versely, if V and W are two decorrelation matrices satisfying 
F = V T V = W T W, then V = OW is an orthogonal rotation 
of W, since (W^-^VW -1 = (VW- 1 ) T VW- 1 = 1 shows 
that = VW _ is an orthogonal matrix. 
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3.2 Prewhitened power spectrum 

The prewhitened nonlinear power spectrum X(k) is the 
Fourier transform of the prewhitened nonlinear correlation 
function X(r) defined in terms of the nonlinear correlation 
function £(r) by (Paper 3 §5) 



X(r) 



2 CM 



(6) 



In this paper, as in Paper 3, hats are used to denote estima- 
tors, so that X(k) with a hat on denotes an estimate of the 
prewhitened power spectrum measured from a galaxy survey 
(Paper 3 §7). The quantity X(k) without a hat denotes the 
prior prewhitened power spectrum, the model power spec- 
trum whose viability is tested in a likelihood analysis. 

The covariance of the prewhitened nonlinear power 
spectrum is approximately equal to the inverse of its Fisher 
matrix, denoted E alS (Paper 3 §6), 

= E af3 . 



(AX a AXp 



(7) 



The weighting of data that yields the best (minimum 
variance) estimate of power depends on the choice of prior 
power, so that the estimated prewhitened power X(k) de- 
pends (weakly) on the prior prewhitened power X(k). An 
incorrect guess for the prior will not bias the estimator high 
or low, but merely makes it slightly noisier than the mini- 
mum allowed by the Fisher matrix (Tegmark et al. 1997). 
The Maximum Likelihood solution to the prewhitened power 
spectrum of a galaxy survey can be obtained by folding the 
estimate X(k) back into the prior X(k) and iterating to 
convergence (Tegmark et al. 1998). 

In this paper, the prior power X(k) is simply treated 
as some fiducial prewhitened power spectrum. In all exam- 
ples illustrated, the prior power spectrum is taken to be an 
observationally concordant ACDM model from Eisenstein & 
Hu (1998). 

3.3 Band-powers 

It is desirable to define estimated band-powers B(k a ) 
and band-power windows W{k a ,k0) for the estimated 
prewhitened power spectrum X(k) in a physically sensible 
way. One of the problems is that the prewhitened power 
spectrum is liable to vary by orders of magnitude, so that, 
carelessly defined, a windowed power could be dominated 
by power leaking in from the wings of the band, rather than 
being a sensible average of power in the band. 
Let band-powers B(k a ) be defined by 



B(k a ] 



W(k a ,kg 



X{ky) inkjdkfs 



(8) 



where xW is some scaling function, to be chosen below 
in equation ([l3|). In manipulations, it can be convenient to 
treat \ as a matrix that is diagonal in Fourier space with 
diagonal entries x(k)- The band-power windows W(k a ,k/3) 
in equation (H) are normalized to unit integral 



W(k a ,kp) 4 ^ k ^ 3/2 dlnfcg = 1 , 



(9) 



the slightly complicated form of which is chosen so as to 
make its discretized equivalent, equation ([n]), look sim- 
ple. Continuous matrices must be discretized in order to 



manipulate them numerically. As described in §2.3 of Pa- 
per 3, discretization should be done in such a way as to 
preserve the inner product in Hilbert space. Discretized 
on a logarithmic grid of wavenumbers, a continuous vector 
such as X(k a ) becomes the discrete vector Xk a = X(k a ) 
[47rfc 3 A In fc/(27r) 3 ] 1,/2 , while a continuous matrix such as 



W(k a , kp) becomes the discrete matrix Wfc a 



W{k a ,k ) 



47r(fc a fc /3 ) 3/2 Alnfc/(27r) 3 . Thus, discretized on a logarith- 
mic grid of wavenumbers, equation (^) translates into 
an equation for the discrete band-powers B ka = B(k a ) 
[47rfe 3 Alnfc/(27r) 3 ] 1/2 , 



5>* 



(10) 



Each row of the discrete matrix Wk^kp represents a band- 
power window for the discrete band-power B ka , and it is 
these discrete band-power windows \Nk a k« that are plotted 



in Figures 2M6l Normalizing each discrete band-power win- 



dow Wfe Q fc^ to unit sum 



= i 



(ii) 



ensures that the band-power B ka represents an average of 
prewhitened power in the band. 

The scaling function x(^) m equations (^) or ( |Io| ) is 
introduced to ensure a physically sensible definition of the 
band-power windows. Rewriting equation (^) as 



B, = 



]x(k a )W kaks x(kf3) 1 X k 



(12) 



makes it plain that choosing different scaling functions \ 
is equivalent to rescaling the band-power windows as W — > 

xw*- 1 . 

Suppose for example that the scaling function in equa- 
tion (hol) were chosen to be one, x(^) — 1> so that B = WX. 
Then a band-power Bfe Q could be dominated by power leak- 
ing in from wavenumbers kp where Xfe^ is large, even though 
the window YJk^kp were small there, which is not good. 

A better choice, the one adopted in this paper, is to 
set the scaling function x(k) equal to the discretized prior 
prewhitened power Xk, 



X(k) = X fc , 

so that the discrete band-powers Bfe Q 
defined by 



(13) 



equation (p^|), 



Bk« 
Xfc Q 



w fc 



Xa 



(14) 



Since the expectation is that Xfe^/X^^ « 1, the defini- 
tion ([y]) ensures that the contribution from power X^ at 
wavenumber kp to a band-power Bk a is large where the win- 
dow Wfc^fc^ is large, and small where the window is small, 
as is desirable. 



3.4 Decorrelated band-powers 

Equation ( |l4| ) expresses the estimated band-powers B as lin- 
ear combinations of scaled prewhitened powers X/X. The 
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discrete Fisher matrix G of the scaled prewhitened power 
is (again, it is convenient to treat the scaling function 
as a matrix that is diagonal in Fourier space, with diagonal 
entries Xfc) 

G = X(AXAX T )~ 1 X = XEX (15) 

where Efc afc/3 = E(k a , kp) 47r(fc a fc )3 ) 3/ ' 2 Aln fc/(27r) 3 is the 
discretized Fisher matrix of the prewhitened power (Paper 3 
§6). 

All the band-power windows W constructed in this pa- 
per are decorrelation matrices, satisfying 

G = W T AW (16) 

where A is diagonal in Fourier space. By construction, the es- 
timated scaled band-powers B/X = WX/X are uncorrelated, 
i.e. their covariance matrix is diagonal in Fourier space 

X _1 (ABAB T )X -1 = A _1 . (17) 

The Fisher matrix of the decorrelated scaled band-powers 
B/X is the diagonal matrix 

X(AB AB T } _1 X = A . (18) 

The decorrelation process decorrelates not only the 
scaled prewhitened power X/X, but also the prewhitened 
power spectrum X itself, as is to be expected since the scal- 
ing factors 1/X are just constants. While the band-power 
windows for the scaled prewhitened power are W, the band- 
power windows for the prewhitened power itself are XWX -1 , 
since B = (XWX" 1 )^. The band-power windows XWX -1 
are themselves decorrelation matrices, equation for the 
prewhitened power X, satisfying 

E = (XWX" 1 ) T (X _1 AX _1 )(XWX _1 ) . (19) 

The decorrelated band-powers B = (XWX _1 )X, are uncorre- 
lated because their covariance matrix is diagonal in Fourier 
space (recall again that the scaling function X is effectively 
a diagonal matrix in Fourier space) 

(ABAB T ) = (XWX _1 )(AXAX T )(XWX _1 ) T = XA _1 X.(20) 

The Fisher matrix of the decorrelated band-powers B is the 
diagonal matrix 

(ABAB T ) _1 = X _1 AX _1 . (21) 




.001 .01 .1 1 10 100 1000 



Wavenumber k (/iMpc ) 

Figure 1. Fisher information matrix G of the scaled prewhitened 
power spectra X/X of (top) SDSS and (bottom) LCRS. The Fisher 
matrix was computed in the FKP approximation, in the manner 
described in §6 of Paper 3. The vertical arrow in the lower panel 
indicates the wavenumber k = it/7.5 TiMpe , where the thick- 
ness of a 1?5 slice at median depth equals a half- wavelength, below 
which the FKP approximation may fail in LCRS. The assumed 
prior power spectrum is an observationally concordant ACDM 
model from Eisenstein & Hu (1998), nonlinearly evolved accord- 
ing to the formula of Peacock & Dodds (1996). Each line rep- 
resents a row (or column, since the matrix is symmetric) of the 
Fisher matrix; each row peaks at, or at small k near, the diag- 
onal. Lines are dashed where the Fisher matrix is negative. The 
resolution is A log A: = 1/128. 



3.5 Interpretation of scaled powers as log powers 

It is interesting, although peripheral to the central thread of 
this paper, to note that the scaled powers B/X and X/X can 
be interpreted in terms of log powers, at least in the limit 
of a large quantity of data, where AX eX-X<X. 

In this limit (B/X) - 1 « ln(B/X) and (X/X) - 1 » 
In (X/X). Thus equation @, with one subtracted from both 
sides, can be rewritten 




or equivalently 

In B ka = C ka + ^ Wfe afc/3 In X fc(3 (23) 



where Ck a = Wfc o fc 3 ln(Xfc Q /X^ ) are constants that 

are zero if the prior prewhitened power Xfc happens to vary 
as a power law with wavenumber k, and in practice should be 
close to zero as long as the band-power windows are narrow. 

Equation (|2^) shows that, modulo the small constant 
offsets C, the log band-powers In B can be regarded as win- 
dowed averages of the log prewhitened powers In X. Irrespec- 
tive of the offsets C, the Fisher matrix of the log prewhitened 
powers approximates the Fisher matrix G, equation jlB]), 

(AlnX AlnX T ) _1 X(AXAX T ) _1 X = G (24) 

and the decorrelation matrices W, equation (|l^), can be 
regarded as decorrelation matrices for the log prewhitened 
power lnX. 
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3.6 Fisher matrices of power in SDSS and LCRS 

Constructing decorrelated band-powers requires knowing 
the Fisher matrix. Figure ^ shows the discrete Fisher ma- 
trices G, equation (^), of the scaled prewhitened nonlinear 
power spectra X/X of SDSS and LCRS. The Fisher matrix 
was computed in the FKP approximation, as described in 
§6 of Paper 3. 

The prior power spectrum X(k) is taken to be a ACDM 
model of Eisenstein & Hu (1998), with observationally con- 
cordant parameters Q,a = 0.7, Q. m = 0.3, Qth 2 = 0.02, and 
h = _Ho/(100 kms _1 Mpc _1 ) = 0.65, nonlinearly evolved ac- 
cording to the Peacock & Dodds (1996) formula. iV-body 
simulations by Meiksin, White & Peacock (1999) indicate 
that nonlinear evolution tends to suppress baryonic wiggles 
in the (unprewhitened) power spectrum, whereas the Pea- 
cock & Dodds transformation preserves the wiggles. For sim- 
plicity we retain the Peacock & Dodds formalism, notwith- 
standing its possible defects. 

The FKP approximation is valid at wavelengths small 
compared to the scale of the survey, so should work better at 
smaller scales in surveys with broad contiguous sky coverage, 
and worse at larger scales in slice or pencil beam surveys. 
Thus the FKP approximation should work reasonably well 
at all but the largest scales in SDSS, but probably fails in 
LCRS at intermediate and large scales. We suspect that the 
FKP approximation is liable to underestimate the Fisher in- 
formation (overestimate the error bars) in LCRS, since the 
density in a thin slice is correlated with (and hence contains 
information about) the density outside the survey volume. A 
reliable assessment of the extent to which the FKP approx- 
imation under- or over-estimates information awaits future 
explicit calculations, such as described in Tegmark et al. 
(1998). 

At the median depth at 300/i _1 Mpc of the LCRS sur- 
vey, each 1?5 slice is 7.5fa _1 Mpc thick, so the FKP ap- 
proximation might be expected to break down in LCRS at 
wavenumbers k ty/7.5 /iMpc -1 = 0.4/iMpc -1 . 

The Fisher matrix shown in Figure [l contains nega- 
tive elements, notably at intermediate wavenumbers k ~ 
l/iMpc -1 , whereas for Gaussian fluctuations the Fisher 
matrix would be everywhere positive. According to equa- 
tion (32) of Paper 3, the Fisher matrix of the power spec- 
trum for Gaussian fluctuations would be 



(25) 



F oc/3 = >_ D * c -Uk c -ljl D ^ _ 

Thus in Fourier space the elements of the Fisher matrix of 
the power spectrum 



F(k a ,kp) 



\C (k a rii, hprij)] 



(4tt) 2 



(26) 



would all be necessarily positive. The quantities dfi in equa- 
tion ( pr^ ) are intervals of solid angle, and the integration is 
over all unit directions n of the wavevectors. This positivity 
of the Fisher matrix for Gaussian fluctuations is separate 
from and in addition to the fact that the Fisher matrix is 
positive definite, i.e. has all positive eigenvalues. 

The fact that some elements of the Fisher matrix in Fig- 
ure |l| are negative is a consequence of nonlinearity. Nonlinear 
evolution induces broad positive correlations in the covari- 
ance matrix of estimates of the power spectrum (Paper 3 
Fig. 2), leading to negative elements in the Fisher matrix, 
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Figure 2. Principal component decomposition of the Fisher ma- 
trix G of the scaled prewhitened power spectrum X/X of SDSS, 
for a ACDM prior power spectrum. Each band-power window is 
an orthonormal eigenfunction of G, or equivalently an eigenfunc- 
tion of the covariance X _1 (AX AX T )X" 1 . The eigenfunctions are 
labelled in order of the expected variances of the corresponding 
band-power, the lowest numbered eigenfunction having the small- 
est variance. The Fisher matrix here was computed over the range 
k = 0.01-100 /iMpc -1 , at a resolution of Alogfc = 1/32. 



the inverse of the covariance of power. Prewhitening the non- 
linear power spectrum substantially narrows the covariance 
of power, and diminishes the amount of negativity in the 
Fisher matrix, but negative elements remain. 



4 EXAMPLE BAND-POWERS 

4.1 Principal Component Decomposition 

Let be an orthogonal matrix that diagonalizes the Fisher 
matrix G of the scaled prewhitened power X/X 



G = 0'AO 



(27) 



Then is a decorrelation matrix. The decorrelated band 
powers OX/X constitute the principal component decompo- 
sition of the scaled prewhitened power spectrum. 

Figure ^| shows the band-power windows, the rows of 
O, for the principal component decomposition of the scaled 
prewhitened power spectrum of SDSS. Since the Fisher ma- 
trix G is already fairly narrow in Fourier space, Figure |l|, one 
might have thought that its eigenmodes would be similarly 
narrow in Fourier space. But this is not so. In fact the eigen- 
modes, the band-power windows, are generally broad and, 
aside from the first, fundamental mode, generally wiggly. 
This makes the principal component decomposition of the 
power spectrum of little practical use, as previously con- 
cluded in Paper 2. This should of course not be confused 
with a principal component decomposition of the density 
field itself, which can be of great utility (Vogeley & Szalay 
1996; Tegmark et al. 1998). 

Elsewhere in this paper plotted band-power windows 
are scaled to unit sum over the window, equation (|ll"|), but 
in Figure |i| the band-power windows are scaled instead to 
unit sum of squares 
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E 



(28) 



In other words, the band-power windows plotted in Figure ^ 
are just the orthonormal eigenfunctions of the Fisher matrix 
G. The problem with scaling the band-power windows to 
unit sum is that it makes the wigglier windows oscillate more 
fiercely, obscuring the less wiggly windows. However, the 
labelling of the band-power windows in Figure ^ is in order 
of the expected variances of the band-powers normalized to 
unit sum over the window, Ofe a fc /3 = 1, not to unit sum 
of squares. 



4.2 Cholesky Decomposition 

Clearly it is desirable to use the infinite freedom of choice 
(§3.1) of decorrelation matrices to engineer band-power win- 
dows that are narrow and not wiggly. To understand why 
principal component decomposition works badly, consider 
the mathematical theorem that the eigenmodes of a real 
symmetric matrix are unique up to arbitrary orthogonal ro- 
tations amongst degenerate eigenmodes, that is, amongst 
eigenmodes with the same eigenvalue. In the present case 
the eigenvalues of the Fisher matrix, even if not degenerate, 
are nevertheless many and finely spaced, and there is much 
potential for nearly degenerate eigenmodes to mix in random 
unpleasant ways. This suggests that mixing can be reduced 
in Fourier space by lifting the degeneracy of eigenvalues in 
Fourier space. 

One way to lift the degeneracy is to multiply the rows 
and columns of the Fisher matrix by a strongly varying func- 
tion 7(fe) of wavenumber k before diagonalizing. That is, 
consider diagonalizing the scaled Fisher matrix 7 T G7 



7 T G 7 = T A0 



(29) 



where 7 is some scaling matrix, is orthogonal, and A is 
diagonal. Then W = 07" 1 is a decorrelation matrix, with 
G = W T AW. Now take the scaling matrix 7 to be a diagonal 
matrix in Fourier space with diagonal entries y(k), and let 
7(fe) be strongly varying with k. 

In the limit of an infinitely steeply varying scaling func- 
tion, j(k{) ^> 7(^2) 3> • • ., the resulting decorrelation ma- 
trix W becomes upper triangular, as argued in §5 of Pa- 
per 2. As pointed out by Tegmark & Hamilton (1998) and 
Tegmark (1997a), this choice of decorrelation matrix is just 
a Cholesky decomposition of the Fisher matrix 



U 1 U 



(30) 



where U is upper triangular. 

If the scaling function 7 is chosen to be infinitely steep 
in the opposite direction, 7(fa) <C 7 (fe) <C • • ., then the 
resulting decorrelation matrix is lower triangular. This is 
equivalent to a Cholesky decomposition of the Fisher matrix 



L L 



(31) 



with L lower triangular. 

Such a Cholesky decomposition has been successfully 
employed to construct a decorrelated power spectrum of the 
CMB from both the COBE DMR data (Tegmark & Hamil- 
ton 1998) and the 3 year Saskatoon data (Knox, Bond & 
Jaffe 1998). 

Figure |^ shows band-power windows from the upper 
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Figure 3. Sample of band-power windows from the Cholesky 
decomposition of the Fisher matrix G of the scaled prewhitened 
power spectrum of (top) SDSS and (bottom) LCRS. The Fisher 
matrix here was truncated to a maximum wavenumber of fc max = 
0.015 /iMpc" 1 in SDSS, and fe max = 0.034 h Mpc" 1 in LCRS, 
to ensure that all its eigenvalues were positive. The resolution 
is Alogfc = 1/32, and only every eighth band-power window is 
plotted. The band-power window at the smallest wavenumber is 
also shown multiplied by a factor of 10, as a dotted line. The 
vertical arrow on the LCRS plot indicates the wavenumber k = 
7r/7.5 /iMpc -1 below which the FKP approximation may fail in 
LCRS. 



triangular Cholesky decomposition of the Fisher matrices G 
of the scaled prewhitened power of SDSS and LCRS. These 
band-powers are a considerable improvement over the princi- 
pal component decomposition of Figure in the sense that 
they are narrower and less wiggly. However, the Cholesky 
decomposition has two defects. The first defect is that the 
band-power windows are skewed. By construction, the upper 
triangular Cholesky windows vanish to the left of the diag- 
onal element. Figure |^ shows that the Cholesky windows 
have a tail to the right that, although small, is neverthe- 
less large enough to start to become worrying, especially in 
LCRS. A lower triangular Cholesky decomposition leads to 
band-power windows skewed in the opposite direction (not 
plotted). 

The second defect of the Cholesky decomposition is that 
it does not tolerate negative eigenvalues in the Fisher ma- 
trix well, a point previously discussed in Paper 2. In Fig- 
ure ^, the Fisher matrices have been truncated to a max- 
imum wavenumber of fe max = 0.015 ft Mpc -1 in SDSS, and 

"-max — 

0.034/1 Mpc" 1 in LCRS, to ensure that the com- 
puted Fisher matrix has no negative eigenvalues. Although 
the Fisher matrix should in theory be positive definite, 
meaning that all its eigenvalues should be positive, numeri- 
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cally it acquires negative eigenvalues when the resolution in 
k is increased to the point where the spacing Ak between 
adjacent wavenumbers is less than of the order of the inverse 
scale size of the survey. Presumably negative eigenvalues ap- 
pear thanks to a combination of numerical noise and the 
various approximations that go into evaluating the Fisher 
matrix. The band-power windows associated with negative 
eigenvalues are ill-behaved, unlike those shown in Figure ^| 
The ill-behaved band-power windows tend to infect their 
neighbours, and as more and more eigenvalues become nega- 
tive with increasing resolution, the whole system of Cholesky 
windows devolves into chaos. 

Replacing the Fisher matrix with a version of it in which 
negative eigenvalues are replaced by zero does not solve the 
difficulty. 

It is possible to extend the band-powers to smaller 
wavenumbers by using a linear instead of logarithmic bin- 
ning in wavenumber. With linear binning, the maximum 
wavenumber is fc max = 7r/1024 /iMpc -1 in SDSS, fc max = 
tt/640 /iMpc" 1 in LCRS. 



4.3 Square Root of the Fisher Matrix 

While the Cholesky band-powers are a definite improvement 
over the principal component decomposition, evidently one 
would like the band-power windows to be more symmetric 
about their peaks. One way to achieve this is to choose the 
decorrelation matrix to be symmetric, which corresponds to 
choosing it to be the square root of the Fisher matrix 



G = G 1/Z G 



l/2 r l/2 



(32) 



The band-powers G 1 ^ 2 X/X are uncorrelated, with unit vari- 
ance. Tegmark & Hamilton (1998) previously used the 
square root of the Fisher matrix to construct the decorre- 
lated power spectrum of CMB fluctuations from the COBE 
data, and found the resulting window functions to be nar- 
row, non-negative and approximately symmetric. One might 
therefore hope that these properties would hold also for the 
case of galaxy surveys. 

Figure |i| shows the band-power windows constructed 
from the square root of the Fisher matrix of the scaled 
prewhitened power of SDSS and LCRS. The plotted win- 
dows, the rows of G 1 ^ 2 , are normalized to unit sum as 
in equation (|ll|). The band-power windows are all nicely 
narrow and symmetrical. The Fisher matrix G itself is al- 
ready fairly narrow about the diagonal, Figure [l], and tak- 
ing its square root narrows it even further, in much the 
same way that taking the square root of a Gaussian matrix 



-(x-x') 2 



narrows it by 2 1//2 , to M 



1/2 



-2(x-x') 2 



In practice, the square root of the Fisher matrix is con- 
structed by diagonalizing the Fisher matrix, G = T A0, 
and setting 



-1/2 



= T A 1/2 



(33) 



the positive square root of the eigenvalues b eing taken. Neg- 
ative eigenvalues, which as remarked in £ 4.2 presumably 



arise from a combination of numerical noise and the approx- 
imations that go into evaluatin g th e Fisher matrix, must be 
replaced by zero in equation (B3[). The negative eigenval- 
ues are invariably small in absolute value compared to the 
largest positive eigenvalue. 




.01 .1 1 

Wavenumber k (/iMpc -1 ) 

Figure 4. Sample of band-power windows from the square root 
G 1 / 2 of the Fisher matrix of the scaled prewhitened power spec- 
trum of (top) SDSS and (bottom) LCRS. The range of the ma- 
trix is k = 0.001-100 hMpc -1 , with resolution Alogfc = 1/32. 
Every eighth band-power window is plotted, starting at k = 
0.0042 /iMpc" 1 in SDSS, and k = 0.0066 h Mpc -1 in LCRS. For 
SDSS, the band-power window at the smallest wavenumber, a 
nominal wavenumber of k = 0.001 /iMpc -1 , is shown multiplied 
by a factor of 10, as a dotted line. What looks like a shaded curve 
(but is actually just a finely oscillating line) underlying this last 
band-power is effectively the same band-power, at the same nomi- 
nal wavenumber of k = 0.001 /iMpc -1 , but computed at 32 times 
higher resolution, Alogfc = 1/1024, over the range k = 0.001— 
0.01 /iMpc - 1 (the range was restricted to keep the computational 
size down). The high-resolution band-power is multiplied by 320 
to raise it to the same level as the low-resolution band-power (re- 
call that the band-powers here are normalized so that the sum 
of their elements is one). The vertical arrow on the LCRS plot 
indicates the wavenumber k = tt/7.5 /iMpc" 1 below which the 
FKP approximation may fail in LCRS. 



Replacing negative eigenvalues in equation (tS3|) by ze- 
ros causes the band-power windows, the rows of G~, to be- 
come linearly dependent, and the resulting band-powers to 
become correlated. Remarkably enough, however, it is still 
possible to regard the band-powers G 1//2 X/X as remaining 
uncorrelated with unit variance. The negative eigenvalues 
correspond to noisy eigenmodes of the Fisher matrix. Re- 
placing the negative eigenvalues by zero in equation ( p3| ) 
corresponds to eliminating the contribution of these noisy 
modes from the band-power windows G 1 ^ 2 . Now if the Fisher 
matrix were calculated perfectly, then it would have no neg- 
ative eigenvalues, so correctly the noisy modes would make 
some positive contribution to G 1 . The contribution would 
however be small because the modes are noisy and the cor- 
responding eigenvalues are small. Thus one should imagine 
that the G 1//2 that results from setting negative eigenvalues 
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to zero differs oniy siightiy from the true G 1 ^ 2 that yields 
band powers G 1//2 X/X with unit covariance matrix. 

As the resolution of the Fisher matrix is increased, more 
and more eigenvalues become negative. The effect on the 
band-power windows is intriguing and enlightening. If for 
example the resolution is doubled, then there are twice as 
many band-powers. The old band-power windows retain es- 
sentially the same shapes as before, except that they are 
computed with twice as much resolution. The new band- 
power windows interleave with the old ones, and have shapes 
that vary smoothly between the shapes of the adjacent band- 
powers. Normalized to unit sum over the window, the vari- 
ance of each band-power doubles, in just such a fashion that 
the net information content, the summed inverse variance, 
remains constant. This behaviour persists to the highest res- 
olution that we have tested it, Alogfc = 1/1024. 

It is worth emphasizing just how remarkable this be- 
haviour of the band-powers is. It seems as though the decor- 
related band-power windows are converging, as the reso- 
lution goes to infinity, towards a well-defined continuous 
shape. Yet, in spite of the asymptotic constancy of shape, 
the variance associated with each window seems to continue 
increasing inversely with resolution, i.e. inversely with the 
number of windows, tending to infinity in the continuous 
limit in such a fashion that the net variance over any fixed 
interval of wavenumber remains constant. 

How can this be? How can there be infinitely many 
uncorrelated band-power estimates? It would seem that as 
the resolution increases, the band-power windows manage to 
remain decorrelated by incorporating into themselves small 
contributions of noisy modes that increase the variance of 
the band-power while changing the windows only slightly. 
In the continuum limit, the variance of each band-power 
would continue to increase indefinitely while the shape of 
the window changes only infinitesimally. 

In practice, the band-power windows at the largest 
scales, nominal wavenumbers less than the natural resolu- 
tion of the survey, become jittery at high resolution. The 
'natural resolution' here is defined empirically, as the high- 
est resolution for which all the eigenvalues of the Fisher ma- 
trix remain numerically positive, Afc = 7r/1024 /iMpc -1 in 
SDSS, and Afc = tt/640 /iMpc" 1 in LCRS. We attribute the 
jitter in part to the fact that the pair integrals R(r; fi) used 
to compute the Fisher matrix are themselves only measured 
with finite resolution and accuracy, and in part to the fact 
that there is practically no signal at the largest scales, so all 
the relevant eigenvalues are small, and the numerics have a 
harder time distinguishing signal from noise. Perhaps there 
is a better algorithm than the simple one adopted here of set- 
ting negative eigenvalues to zero, but the simple algorithm 
does seem to work well enough. 

The situation is illustrated in Figure ||, which shows 
the first band-power window in the sequence for SDSS, 
the one nominally corresponding to k = 0.001 h Mpc , 
multiplied by a factor of 10 to show it more clearly. This 
nominal wavenumber is smaller by a factor of ~ 3 than 
the smallest measurable wavenumber fc ~ 0.003 ft Mpc -1 in 
SDSS, and the computed band-power accordingly retreats 
to a larger effective wavenumber, where there is signal. At 
scales comparable to the natural resolution of the survey, 
k ~ Afc ~ 0.003/1 Mpc" 1 , the resolution Alogfc = 1/32 
used in Figure ^ is 10 times higher than the natural resolu- 
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Figure 5. Sample of band-power windows from the square root 
G 1 / 2 of the Fisher matrix of the scaled prewhitened power spec- 
trum of (top) SDSS and (bottom) LCRS, with linear rather than 
logarithmic gridding in wavenumber. The resolution is Afc = 
7r/4096 /iMpc -1 , finer than the natural resolution of the sur- 
vey by a factor of about 4 in SDSS, 6 in LCRS. Every 32nd 
band-power is plotted, starting at a nominal wavenumber of 
k = tt/1024 h Mpc -1 . For LCRS, the first band-power plotted, at 
nominal wavenumber k = 7r/1024 ftMpc - , is jittery. For LRCS, 
the FKP approximation may fail at all wavenumbers shown. 



tion. For comparison, Figure ^] also shows (shaded) the same 
band-power computed at a resolution of Alogfc = 1/1024, 
several hundred times the natural resolution. At this high 
resolution, the band-power oscillates finely, but overall re- 
mains under control. Such robust behaviour contrasts with 
the chaotic behaviour of the Cholesky windows when simi- 
larly pushed beyond the natural resolution of the survey. 

The methods of Paper 3 permit the Fisher matrix to 
be discretized over any arbitrary grid of wavenumbers (al- 
though all the explicit examples in Paper 3 employ a loga- 
rithmic grid) . Figure |E] shows the band-power windows con- 
structed from the square root of the Fisher matrix of the 
scaled prewhitened power of SDSS for a linearly spaced 
rather than logarithmically spaced grid of wavenumbers. 
These band-power windows have essentially the same shapes 
as those computed on the logarithmic grid, Figure^. The dif- 
ference is in the number of band-power windows and in the 
resolution with which they are defined. 

Figure ^| shows, more clearly than Figure ^, that the 
band-powers are broader for LCRS than for SDSS, reflect- 
ing the smaller effective volume, hence lower resolution in 
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Figure 6. Comparison of band-power windows obtained from 
(thin line) the square root of the Fisher matrix, G 1 / 2 , versus (thick 
line) the scaled square root of the Fisher matrix, (7 T G 7) 1 / 2 7~ 1 
with scaling function "/(k) = k 3 ' 2 , for the case of a perfect, noise- 
less survey. Only the single band-power centred at k = IMpc -1 
is shown. Lines are dashed where the windows are negative. The 
resolution is Alogfe = 1/128. The figure illustrates that in cases 
where the Fisher matrix varies steeply, scaling the Fisher matrix 
before taking its square root can help to produce more symmetric 
band-powers. 



Figure 7. Inverse variance, or information, per logarithmic in- 
terval of wavenumber in the scaled prewhitened power spectra 
of the SDSS, LCRS, and IRAS 1.2 Jy surveys, for the ACDM 
prior power spectrum. This information is the inverse variance 
on each band-power, if the band-powers are an e-fold wide. The 
information in the LCRS is probably underestimated here, espe- 
cially at larger scales, because of the assumption of the FKP 
approximation. The vertical arrow indicates the wavenumber 
k = 7r/7.5 hMpc -1 below which the FKP approximation may 
fail in LCRS. 



Fourier space, of LCRS. The slice geometry of LCRS leads 
to wings on the band-powers. The wings are fairly mild 
here, but in pencil-beam surveys the wings can become quite 
broad, potentially leading to significant aliasing between 
power at small and large wavenumbers (Kaiser & Peacock 
1991). 

For LCRS, the first band-power plotted in Figure H, at 
a nominal wavenumber of k = 7r/1024 ftMpc" 1 , is jittery, 
illustrating again that band-powers at nominal wavenum- 
bers less than the natural resolution of the survey, k — 
7t/640 /iMpc -1 for LCRS, tend to become jittery at high 
resolution. Again, we attribute this in part to the fact that 
the pair integrals R(r; /i) used to compute the Fisher matrix 
are themselves computed with finite accuracy (in LCRS, the 
pair integrals were evaluated by the Monte Carlo method), 
and in part to the fact that the numerics have a harder time 
distinguishing signal from noise when the signal is weak. 



4.4 Scaled Square Root of the Fisher Matrix 

Notwithstanding the success of the square root of the Fisher 
matrix in many cases, it does not work perfectly in all cases. 
The problem is that the symmetry of G 1 ^ 2 does not imply 
symmetry of the band-power windows, the rows of G 1//2 , 
about their diagonals. If the Fisher matrix G varies steeply 
along the diagonal, then the band-power windows can turn 
out quite skewed. In practice, this happens for example in 
the case of a perfect, shot-noiseless survey, where the infor- 
mation contained in the power spectrum increases without 
limit as k — > oo. The case of a noiseless survey was applied 
in §9.3 of Paper 3 to compute the effective FKP constants 
/i(fc) to be used in an FKP pair-weighting when measuring 
the prewhitened nonlinear power spectrum of a survey. 



The difficulty can be remedied by scaling the Fisher 
matrix before taking its square root. Let 7 be any scaling 
matrix. Then 



W= ( 7 T G 7 ) 1/ V 1 



(34) 



is a decorrelation matrix, satisfying G = W T W. 

In the case of the perfect, noiseless survey, a diagonal 
scaling matrix 7 with diagonal elements 



7 (fc) = fc 3/2 



(35) 



proves empirically to work well. Figure [| compares band- 
power windows obtained from the scaled versus unsealed 
square root of the Fisher matrix, for the perfect, noiseless 
survey. For clarity only a single band-power window, the 
one centred at k = 1 /iMpc -1 , is shown, but this window is 
representative. Scaling the Fisher matrix with j(k) = k 3 ^ 2 
before taking its square root in this case helps to rectify the 
skew in the band-power windows from the unsealed square 
root G 1//2 . The results shown in Figure 11 of Paper 3 were ob- 
tained using the decorrelation matrix of equation (|34|) scaled 
by the scaling function 'y(k) = k 3 ^ 2 , equation (pS5|). 



5 DECORRELATED POWER SPECTRUM 

The goal of this paper has been to show how to obtain uncor- 
rected band-po wers B that can be plotted, with error bars, 
on a graph. In §4.3 it was found that, if the band-powers 



are decorrelated with the square root of the Fisher matrix, 
then the band-powers can be computed at as high (or low) a 
resolution as one cares, and the band-powers will remain ef- 
fectively uncorrelated. The higher the resolution, the larger 
the error bars on the band-powers. The same result holds 
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Figure 8. Dccorrelated power spectrum expected from SDSS. 
The prior power spectrum is the observationally concordant 
ACDM model of Eisenstein & Hu (1998) used throughout this 
paper, with parameters as indicated, nonlinearly evolved accord- 
ing to the formula of Peacock & Dodds (1996). Both the lin- 
ear power spectrum §l(&) and the prewhitened nonlinear power 
spectrum X(k) are shown. Error bars on the prewhitened non- 
linear power spectrum are linearly spaced in wavenumber at 
Ak = tt/1024 hMpc" 1 up to k = 0.04hMpc~ 1 , then logarith- 
mically spaced at Alogfc = 1/32. Eisenstein & Hu also provide a 
'nowiggle' power spectrum that is a smooth fit through the bary- 
onic wiggles in the spectrum. The lower panel shows the ratio 
of the prewhitened nonlinear power spectrum to the prewhitened 
nonlinear nowiggle power spectrum. 



if the band-powers are decorrelated with the scaled square 
root of the Fisher matrix, equation (|3~i|). 

If the band-power windows W are scaled to unit sum, 
equation (p"l|), then the Fisher matrix of the scaled band- 
powers B/X is the diagonal matrix A in G = W T AW, equa- 
tion (^). The quantity that remains invariant with respect 
to resolution (and linear or logarithmic binning) is the in- 
verse variance, also called the information, per unit log 
wavenumber d//d In fc 

dlnfc A In A; [ ' 

where A^ is the diagonal element of A associated with the 
scaled band-power B^/X^, and Alnfc is the resolution of the 
matrix at wavenumber fc. 

The association of a band-power B^ with wavenumber 
fc relies on the band-power window being narrow about that 
wavenumber. In practice the band-power windows have a 
finite width, and the correspondence with wave number is 
not exact, a r eflection of the uncertainty principle. As dis- 
cussed in £4.3, the decorrelated band-power windows would 



remain of finite width even if they were resolved at infinite 
resolution. 

Figure |^ shows the information per unit log wavenum- 
ber d//dlnfc in the scaled prewhitened power spectra of 



SDSS, LCRS, and also the IRAS 1.2 Jy survey, for the 
ACDM prior power spectrum. These information curves are 
similar to those presented by Tegmark (1997b), but are more 
accurate, and valid also in the nonlinear regime. 

It should be cautioned that the information plotted in 
Figure ^ comes from a Fisher matrix computed in the FKP 
approximation, which is correct only for scales small com- 
pared to the scale of the survey. The FKP approximation 
tends to misestimate, probably underestimate, the informa- 
tion contributed by regions near (within a wavelength of) 
survey boundaries, since it assumes that those regions are 
accompanied by more correlated neighbours than is actually 
the case. The problem is most severe in surveys like LCRS, 
where everyone lives near the coast. Thus the information 
plotted in Figure [j] probably underestimates the true infor- 
mation in the LCRS, especially at larger scales. 

For a particular choice of resolution in wavenumber, 
the information per unit log wavenumber translates into an 
uncertainty in the corresponding uncorrelated band-powers. 
The amount of information Alt in a scaled band-power 
Bfc/Xfc of width A In A; is equal to the corresponding diagonal 
element Afc of the Fisher matrix of the scaled band-powers: 



Ah 



dl 



Alnfc = A fc 



(37) 



dlnfc 

The inverse square root of this is the expected error in the 
scaled band-power 



<(AB fc /X fc ) 2 > 1/2 = (AI k ) 



-1/2 



= a; 



-1/2 



(38) 



Figure |^ illustrates an example of the error bars ex- 
pected on the decorrelated prewhitened nonlinear power 
spectrum of SDSS, for the ACDM prior power spectrum. As 
always in likelihood analysis, 'the errors are attached to the 
model, not to the data'. Figure |§| demonstrates that, if bary- 
onic wiggles are present at the expected level, and if nonlin- 
ear evolution leaves at least the first wiggle intact, as sug- 
gested by A^-body simulations (Meiksin, White & Peacock 
1999), then SDSS should be able to recover them. According 
to the prognostications of Eisenstein, Hu & Tegmark (1999), 
the detection of baryonic features in the galaxy power spec- 
trum should assist greatly in the business of inferring cos- 
mological parameters from a combination of CMB and large 
scale structure data. 



6 CONCLUSIONS 

Amongst the infinity of possible ways to resolve the galaxy 
power spectrum into decorrelated band-powers, the square 
root of the Fisher matrix, or a scaled version thereof, offers 
a particularly good choice. The resulting band-powers are 
narrow, approximately symmetric, and well-behaved in the 
presence of noise. 

By contrast, a principal component decomposition 
yields band-powers that are broad and wiggly, which renders 
them of little practical utility. A Cholesky decomposition of 
the Fisher matrix works better than principal component 
decomposition, but not as well as the square root of the 
Fisher matrix. On the good side, Cholesky band-power win- 
dows are narrow and not wiggly; on the bad side, Cholesky 
band-power windows are skewed to one side, and they re- 
spond poorly to the presence of small negative eigenvalues 
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in the Fisher matrix, as can occur because of a combination 
of numerical noise and the various approximations that go 
into computing the Fisher matrix. 

In summary, the square root of the Fisher matrix is a 
useful tool for decorrelating the power spectrum not only of 
CMB fluctuations (Tegmark & Hamilton 1998), but also of 
galaxy redshift surveys. 

We conclude with the caveat that this paper, like the 
preceding Paper 3, has ignored redshift distortions, and 
other complicating factors such as light-to-mass bias. In any 
realistic analysis of real galaxy surveys, such complications 
must be taken into account. 
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